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We construct models of interacting itinerant non-Abelian anyons moving along one-dimensional 
chains. We focus on itinerant Ising (Majorana) and Fibonacci anyons, which are, respectively, related 
to SU(2)2 and SU(2)3 anyons and also, respectively, describe quasiparticles of the Moore- Read and 
Zs-Read-Rezayi fractional quantum Hall states. Following the derivation of the electronic large-U 
effective Hubbard model, we derive effective anyonic t-J models for the low-energy sectors. Solving 
these models by exact diagonalization, we find a fractionalization of the anyons into charge and 
(neutral) anyonic degrees of freedom - a generalization of spin-charge separation of electrons which 
occurs in Luttinger liquids. A detailed description of the excitation spectrum can be performed by 
combining spectra for charge and anyonic sectors. The anyonic sector is the one of a squeezed chain 
of localized interacting anyons, and hence is described by the same conformal field theory (CFT), 
with central charge c = 1/2 for Ising anyons and c = 7/10 or c = 4/5 for Fibonacci anyons with 
antiferromagnetic or ferromagnetic coupling, respectively. The charge sector is the spectrum of a 
chain of hardcore bosons subject to phase shifts which coincide with the momenta of the combined 
anyonic eigenstates, revealing a subtle coupling between charge and anyonic excitations at the 
microscopic level (which we also find to be present in Luttinger liquids), despite the macroscopic 
fractionalization. The combined central charge extracted from the entanglement entropy between 
segments of the chain is shown to be 1 + c, where c is the central charge of the underlying CFT of 
the localized anyon (squeezed) chain. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.Mg 



I. INTRODUCTION 

One of the most significant pursuits in condensed mat- 
ter physics is the search for quasiparticles or excitations 
that obey non-Abelian exchange statistics^ The most 
prominent candidates (at present) are excitations of non- 
Abelian quantum Hall states^^. In particular, there is 
evidence from tunneling^^ and interferometry^^'^^ exper- 
iments supporting the existence of non-Abelian quasipar- 
ticle excitations for the v = b/2 quantum Hall state^^^^. 

The leading candidate quantum Hall states to describe 
the electronic ground state of the quantum Hall plateau 
at filling fraction v = b/2 are the Moore-Read (MR) 
Pfaffian state^ or its particle- hole conjugate, the "anti- 
Pfaffian" (aPf) state^'^. Despite the differences between 
the MR and aPf states (which manifest in the detailed 
structure of the edge states), the non-Abelian anyonic 
structure of their bulk quasiparticles are simply the com- 
plex conjugates of each other. Both can be described in 
terms of an Ising-type anyon model. 

One of the leading candidates to describe the experi- 
mentally observed v = 12/b quantum Hall plateau^^'^*" is 
the /c = 3 Read-Rezayi (RR) state^ (a generalization of 
the MR state), or, more precisely, its particle-hole con- 
jugate (RR). The non-Abelian quasiparticles of the RR 
and RR states are of Fibonacci type, and the v = 12/5 
quantum Hall state is the leading candidate system host- 



ing such non-Abelian anyons. The other leading candi- 
date for describing the u = 12/5 quantum Hall effect is 
provided by Bonderson-Slingerland (BS) states^ obtained 
hierarchically from the MR and aPf u = 5/2 states (by 
condensing Laughlin-type quasiholes). The quasiparti- 
cles of these BS states have a similar Ising-type non- 
Abelian structure as their MR and aPf parent states. 
Numerical studies of the u = 12/5 quantum Hall state 
found the RR and BS candidates to be in close competi- 
tion^^. 

Interestingly, a different hierarchical construction over 
the MR state (condensing fundamental non-Abelian 
quasielectrons) produces a candidate state for filling u = 
18/7 that possesses non-Abelian quasiparticles of the Fi- 
bonacci type^^, similar to the ones appearing in the non- 
Abelian spin-singlet (NASS) state^^. However, a quan- 
tum Hall state at u = 18/7 has, so far, not been experi- 
mentally observed. 

Another promising class of candidates for realizing 
non-Abelian quasiparticles is provided by systems with 
the so-called emergent Majorana zero modes, which be- 
have like Ising-type anyons under exchange. Majorana 
zero modes were originally predicted to exist in vortex 
cores of chiral p-wave superconductors^^ '^^ or at the ends 
of one-dimensional polarized superconductors^^. More 
recently, it was shown that Majorana fermions can form 
at the interface of a strong topological insulator and an 
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s-wave superconductor^^. This idea for realizing Majo- 
rana fermion zero modes was further developed by several 
groups^^^^, who proposed similar superconducting het- 
erostructures based on semiconductors exhibiting strong 
spin-orbit coupling, rather than topological insulators. 
For a review, see Ref. 29. Efforts to physically imple- 
ment these latest proposals have been made in recent 
experiments^^^^ consisting of electrons tunneling into a 
nanowire (set-up to be in a topological phase supporting 
edge Majorana modes). However, neither the exponen- 
tial localization to the edges of the observed zero-modes 
nor the non-Abelian exchange statistics has been probed 
yet. 

Inspired by these recent developments and possible re- 
alizations of quasiparticles with non-Abelian statistics, 
we consider the question of what happens if one confines 
mobile non-Abelian quasiparticles to one-dimensional 
(ID) systems. The concept of itinerancy of interact- 
ing non-Abelian quasiparticles is of direct physical sig- 
nificance, and the microscopic models we study can be 
viewed, for example, as (crude) effective models relevant 
to edge modes of quantum Hall and Majorana fermion 
systems. It was established long ago, starting with the 
work of Anderson^^, that electrons confined to one di- 
mension undergo "spin-charge separation," namely the 
electrons falls apart into two pieces, one spinless carrying 
the charge, the other a spinon without charge, carrying 
the spin. These ideas were further developed by several 
people, Tomonaga^^, Luttinger^^, and Haldane^'', who 
introduced the concept of the one-dimensional Luttinger 
liquid. 

In our recent Letter^^, we started to investigate 
the subject of itinerant non-Abelian anyons in a one- 
dimensional system. We established that non-Abelian 
anyons (of which the quasiparticles of the quantum Hall 
states and the Majorana zero modes discussed above are 
prime examples) also undergo a process which resembles 
spin-charge separation. Namely, the non-Abelian anyons 
fractionalize into charge (or density) and anyonic degrees 
of freedom. The model introduced in Ref. 38 was inspired 
by the electronic t-J model^^, which can be viewed as a 
limiting case of the Hubbard model^^^^ - namely in the 
limit of large on-site repulsion - and for which spin-charge 
separation was established analytically at a supersym- 
metric point^^^^ and numerically^^. In this paper, we 
continue our study of itinerant non-Abelian anyons and 
provide greater detail. 

The outline of the paper is as follows: in Sec. II we 
review briefiy general properties of non-Abelian anyons 
in SU(2)/c Chern-Simons theories and, more specifically, 
in non-Abelian quantum Hall states. In Sec. Ill, we show 
that, in close analogy to the electronic case, it is possible 
to: (i) truncate the Hilbert space of the quasiparticles 
of the non-Abelian quantum Hall states confined to a 
one-dimensional geometry (in the case of strong charg- 
ing energy), and (ii) derive low-energy effective anyonic 
t-J models. The charge sectors of the anyonic t-J mod- 
els are derived in Sec. IV. In Sec. V, we present a short 



review of the properties of dense (localized) non-Abelian 
anyon chains. In Sec. VI, we demonstrate that the ex- 
citation spectrum of the anyonic t-J models can be ac- 
curately described by combining spectra for charge and 
anyonic sectors (in a subtle manner), and provide clear 
evidence of the fractionalization of anyons into charge 
and anyonic degrees of freedom. In Sec. VII, we use 
density matrix renormalization group (DMRG) calcula- 
tions to extract the central charge of the Fibonacci t-J 
chain from the entanglement entropy between segments 
of open chains, providing further evidence of fractional- 
ization. In Appendix A, we provide detailed descriptions 
of the related anyon models of the MR and RR states. In 
Appendix B, we describe the details of the quasiparticle 
spectrum truncation for MR and RR anyons used in the 
paper. 



II. NON-ABELIAN ANYONS 
A. General considerations and fusion algebra 

The non-Abelian anyons we will be concerned with in 
this paper can formally be described by SU(2)/e Chern- 
Simons theories, or via a certain quantum deformation 
of SU(2). In either case, the non-Abelian degrees of free- 
doms are captured by the "topological charges" j, which 
can be thought of as "generalized angular momenta." For 
a given SU(2)/c theory, these are constrained to take the 
values j = 0, ^, . . . , |, loosely corresponding to the first 
k representations of SU(2). 

In the same way that the tensor product of SU(2) spins 
can be decomposed into the direct sum of multiplets of 
definite values of J^, one can decompose the product, or 
"fusion" of anyons. This fusion algebra or "fusion rules" 
of a general anyon model takes the form 

axb = Y,KbC (1) 

cec 

where a, 6, and c are topological charge values in the set 
of allowed topological charges C, and the fusion coeffi- 
cients N^^ are non-negative integers indicating the num- 
ber of ways a and b can fuse to produce c. The N^^ must 
be such that the algebra is commutative and associa- 
tive. There must also be a unique "vacuum" or "trivial" 
charge, which we denote as / or 0, for which A^^q = Sac 
The fusion rules of SU(2)/e anyons resemble their SU(2) 
counterpart (but with a finite set of allowed values for 
j and a corresponding truncation of the algebra). In 
particular 

min{ji+j2,fc-ji-j2} 
jl X j2 = ' 

j3 = \jl-j2\ 

where the upper limit is such that the fusion rules are 
associative and obey the constraint that ji < |. 



3 



e {0,1/2,.. .p 



(Xi G {J, cr, i;} 



e e I 

f f Q 



Xi X2 X2 

a a I 
f Q 



Xi X2 X2 

T T I 

m Q Q 



.T| X2 X2 



e 



a 



I 

Q 



(a) 



(b) 



(c) 



FIG. 1: Anyonic fusion trees for (a) electrons, (b) Ising a 
anyons, and (c) Fibonacci anyons. Sites are denoted by the 
(filled or empty) circles at top of the diagrams. Empty cir- 
cles denote vacant sites, which carry the vacuum or trivial 
topological charge or /. The bond labels {xi} encode non- 
local information about the state and their possible values are 
specified for each model. We note that, for Ising anyons, our 
model excludes ip anyons on the sites, but not on the links 
(see the main text). 



We will be particularly interested in the cases k = 2 
and 3, because they are the experimentally most relevant 
non-Abelian anyon models. (The case k = 1 corresponds 
to Abelian anyon models, in which case the Hilbert spaces 
we construct below are one-dimensional, so one can not 
construct non-trivial models.) These cases are related 
to the Ising and Fibonacci anyons models, respectively, 
which we now consider in more detail. 

The anyon model with k = 2 has three anyon types, 
which we will label using the Ising TQFT topological 
charges /, <j, and ip as follows: The vacuum or trivial 
anyon / corresponds to j = 0. The anyon of type a 
corresponds to j = ^. Finally, the (fermionic) anyon 
type i/j corresponds to j = 1. In particular, the fusion 
rules read^^ 



a X a = I -\- t/j a X = a ip x ip = I ^ 



(3) 



in addition to the general relations a x P = P x a and 
I X a = a, which hold in all anyon models, for arbitrary 
a and /3. 

The Fibonacci anyon model correspond to the case k = 
3, where we restrict ourselves to the integer- valued j = 
and 1, which we will label as the vacuum / and the 
Fibonacci anyon r, respectively.^^ The Fibonacci fusion 
rule reads 

r X r = / + r . (4) 

In the following we shall consider itinerant anyons mov- 
ing on one dimensional chains. Pictorial representations 
of such anyonic chains are shown in Fig. l(b,c), together 
with the more familiar case of strongly correlated elec- 
trons shown in Fig. 1(a). 



For the case of strongly correlated electrons, each elec- 
tron carries a unit charge and spin-^. The spin-^ de- 
grees of freedom are taken into account in Fig. 1(a) in 
a slightly unconventional way, utilizing a "fusion tree" 
notation, rather than the usual tensor product of N two- 
dimensional local Hilbert spaces (where N is the number 
of electrons). In this notation, the labels Xi for the links 
of the fusion tree correspond to the total spin obtained 
by combining the spin with that of the ith electron. 
For an open chain, this simply means that Xi is the total 
spin of all the electrons to the left of the label. For a 
periodic chain (in a system on a torus), the label has a 
slightly more abstract interpretation, since the notion of 
all particles to the left or right are not well-defined. We 
use this formulation because it easily generalizes to the 
case of non-Ablian anyons, where there are no local de- 
grees of freedom (i.e. they lack local Hilbert spaces and 
internal quantum numbers, similar to Sz in the case of 
spins). 

The non-Abelian anyons in Fig. l(b,c) may also carry 
electric charge (albeit this typically is a fraction of the 
charge of the electron) , as well as anyonic degrees of free- 
dom. The charge degrees of freedom live on the sites, 
while the bond variables Xi encode the anyonic degrees 
of freedom along the fusion tree, in the same way as the 
labels Xi encoded the spin of the electrons in Fig. 1(a). 
Abelian anyonic degrees of freedom may be treated in the 
same way as electric charge, i.e. locally assigned to the 
sites, since their resulting fusion tree is uniquely deter- 
mined by the local degrees of freedom. The labels Xi are 
not arbitrary, but satisfy the constraint that each triva- 
lent vertex in this fusion tree is permitted by the fusion 
rules. This implies that the size of the internal Hilbert 
space (for a given configuration of particle/anyon posi- 
tions) grows as 2^ in the case of the electrons, (a/2)^ in 
the case of the Ising anyons, and (p^ in the case of the 
Fibonacci anyons, where <p = is the golden ratio. 

Here, N corresponds to the number of electrons, Ising a 
anyons, or Fibonacci anyons. The actual dimension for 
any finite N is, of course, an integer, so these are only the 
leading order scaling (as N 00) for the non-Abelian 
anyons. The sites labeled by / correspond to vacancies, 
and carry no electric charge, spin, or anyonic degrees of 
freedom. 

Before we continue in the next subsection with de- 
scribing the quantum Hall states in which these types of 
anyons are realized, we want to make one remark, which 
will be essential in the subsequent description of the be- 
havior of the itinerant anyons. Despite the fact that we 
will be describing mobile, but identical, anyons, there 
will be a notion of "distinguishability" of the anyons. In 
particular, the various states in the Hilbert space are not 
only characterized by the location of the occupied sites, 
but also by the labels x^, which distinguish the various 
states, given the location of all the anyons. In some sense, 
specifying the precise internal state, corresponding to all 
the anyons as a whole, renders the individual anyons in a 
particular state distinguishable. We will see later on that 
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this seemingly simple observation will play an essential 
role in the effective description of the collective behavior 
of the itinerant anyons. 



B. Non-Abelian quantum Hall states 

We now concentrate on describing the anyonic struc- 
ture of the MR and k = 3 RR states in their fermionic 
incarnations, which are relevant in the electronic quan- 
tum Hall setting. Because of the fermionic nature of the 
states, the anyonic structure is slightly more complicated 
than the 811(2)^ anyons described above. To describe this 
structure, it is best to consider the non-Abelian part sep- 
arately, which is described in terms of Ising anyons for 
the MR state, and Z3 parafermions for the k = 3 RR 
state. 

In the case of the MR state, the non-Abelian sector is 
the Ising theory, whose fields are /, a, ?/^, with the Ising 
fusion rules given above. The quasiparticle types can now 
be specified by the Ising label, together with the electric 
charge. The vacuum is (/, 0), while (a, e/4) is the "funda- 
mental quasiparticle," which, in some sense, carries the 
"smallest" quantum numbers allowed in the MR state, 
i.e. it has the smallest (nonzero) electric charge and re- 
peated fusion generates the entire spectrum of topological 
charges. All other quasiparticles are thus obtained by re- 
peated fusion of this fundamental quasiparticle, using the 
fusion rules above and the additivity of the charge. In ad- 
dition, one needs the rule that quasiparticles which differ 
by fusion of an electron, given by (?/^, e) are to be identi- 
fied. The fact that we identify quasiparticles which "dif- 
fer by a fermion" (or identify the electron with the vac- 
uum) leads to some complications, which are not present 
for the bosonic versions of these quantum Hall states^^, 
but these complications will not concern us here. The re- 
sulting quasiparticle spectrum is given in Fig. 2(b), where 
we have six different quasiparticle types (shown as green 
circles), because (/, e) and (^^,0) are identified, and so 
on. 

In the case of the RR state with /c = 3, the non-Abelian 
structure corresponds to the Za-parafermions, which con- 
ventionally are labeled as /, ?/^2, en, cr2, ^- We will, 
however, use the notation ?/;o = /, tq = £, r2 = (Ji and 
Ti = (72. In this way, the fusion rules take the simple 
form 



Ti X Tj = V^i+j + Ti^j 



(5) 



where the indices are taken modulo 3. This is the direct 
product of the Fibonacci fusion algebra with a Z3 fusion 
algebra. The quasiparticle types can now be specified 
by the Zs-parafermion label, together with the electric 
charge. The vacuum is (V^o^O) while the fundamental 
quasiparticle's quantum numbers are (r2, e/5). As in the 
MR state, the other quasiparticles are obtained by re- 
peated fusion of the fundamental quasiparticle, using the 
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FIG. 2: Charts of the quasiparticle contents of (a) the elec- 
tronic Hubbard model, (b) the MR state and (c) the k — 3 RR 
state. The elementary electric charges of the MR quasiparti- 
cles are multiples of e/4, while that of the RR quasiparticles 
are multiples of e/5. The dark (green) symbols correspond to 
the different particle types. The black dots represents elec- 
trons/holes, which are identified with the vacuum (/, 0) in (b) 
and the vacuum (V^cO) in (c). The grey dots correspond to 
particles which are identified with one of particles correspond- 
ing to a green symbol, as explained in the main text. 



fusion rules above. Quasiparticles which differ by fusion 
of an electron, given by (V^i, e), are to be identified. This 
gives rise to ten different quasiparticle types displayed as 
green circles in Fig. 2(c), with charges 0, . . . ,4e/5 [where 
we note that (ri,e) is identified with (ro,0), and so on]. 

For comparison, we display the relevant quasiparticle 
types in the ordinary electron case in Fig. 2(a): the trivial 
particle (vacant site) /, the electron e, a double occupied 
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FIG. 3: (a) All two-anyon processes can be represented as 
a general "tunneling" term, where topological (and/or elec- 
trical) charge is transferred from one localized quasiparticle 
to the other. Special cases of this include (b) an "interac- 
tion" between the two anyons, for which the localized charges 
are unchanged, and (c) a "hopping" term, where a localized 
anyon moves to a vacant site. 



site d and a spin-1 magnon m. In this case, there is no 
condensate, and hence none of the particles are to be 
identified. 

In the next section, we will describe how we can trun- 
cate the spectrum of particles, in order to come up with 
a tractable model of interacting, itinerant anyons in one 
dimension. 



III. ANYONIC t-J MODELS 
A. Objectives and procedure 

We move now to the construction of the low energy 
models for the itinerant anyons in one dimension, mod- 
eled by a discrete chain. This chain might be a lattice dis- 
cretization of a one-dimensional continuum system such 
as the edge of a quantum Hall liquid or a one-dimensional 
array of quantum dots. On this chain we restrict our- 
selves to short-range (nearest-neighbor) interactions and 
"hopping" terms that can both be expressed as "tunnel- 
ing" processes^'', as sketched in Figs. 3,4. Because the 
anyons are electrically charged, confinement on a quan- 
tum dot or transverse confinement (e.g. in the case of 
edge states) may lead to large Coulombic charging en- 
ergy that strongly discourages multiple occupancy of sites 
and, hence, prohibits the formation of quasiparticle exci- 
tations of larger charge values. 

As seen in the previous section the physical contents of 
the non-Abelian quantum Hall states are very rich and 
we thus want to derive an effective low-energy model, 
similar to the derivation of the t-J model for electrons 
from the Hubbard model. In order to describe the low- 
energy spectra, a well-know strategy consists of build- 
ing up a simpler model by (i) discarding the high-energy 
(quasi)particles and (ii) treating virtual processes involv- 
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FIG. 4: Hopping and interaction terms for (a) electrons, (b) 
Ising anyons, and (c) Fibonacci anyons, expressed in the no- 
tations of Fig. 3. 



ing the fusion of the low-energy (quasi)particles to the 
discarded high-energy (quasi) particles in second-order 
perturbation theory. 



B. Large-U electronic Hubbard model 

To illustrate this procedure, we first take the example 
of a generalized Hubbard model of electrons and show 
how to derive the corresponding t-J model. We start 
from electrons, where the most general SU(2)-symmetric 
single band model with nearest neighbor interactions can 
be written in second quantized notation as 



H 



i ^ 

i 



(6) 



Here ^ and q^ct are the creation an annihilation oper- 
ators of an electron with z-component of spin a, n^^^r = 
c] ^c- ^ are the local spin densities, rii = rii^^ + is the 

total local density, and Si is the spin operator on site i. 

The first term (t) in the Hamiltonian of Eq. (6) is the 
hopping (tunneling) of an electron. The second term 
(Jo) is a spin exchange term, which can be interpreted 
as a two-electron interaction mediated by the tunneling 
of a spin-1 magnon. This term can also be written as 
— JVf~^i^ with Vf~^i being the projector onto the total 
singlet state of two neighboring electrons at sites i and 
i + 1. The third term {V) is a nearest neighbor repulsion, 
which can be interpreted as tunneling of a photon, and 
finally the last term (t/) is the local charging energy. 
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FIG. 5: Schematic energy spectra (arbitrary scale) in the 
presence of a parabohc Coulomb charging energy for the case 
of (a) electrons, (b) Ising anyons, and (c) Fibonacci anyons. 
The Hubbard U energy is shown in (a). Upward arrows show 
the shifts corresponding to the topological contribution to the 
bare energies of the quasiparticles, e.g. A^^ = 1/8 and = 1 
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In this Hubbard model of electrons, we consider three 
different types of "quasiparticles" at the lowest energies: 
the "trivial particle" / (i.e. a vacant site), the electron e, 
and the "doublon" (i, corresponding to a double (elec- 
tronic) occupancy on a site. In Fig. 5(a), we model 
the energy costs for a given quasiparticle type using a 
parabolic Coulombic charging energy of the form 



VcouM) ^ {q- Qc) 



(7) 
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FIG. 6: (a,b) Second-order exchange processes for two elec- 
trons on nearest neighbor sites via the virtual creation of a 
"vacuum" quasiparticle / (i.e. a vacancy) and a "doublon" 
d. These two exchange diagrams can be replaced by a first- 
order magnon exchange diagram, shown in (c), leading to a 
renormalization of the magnon exchange interaction. 



where q is the quasiparticle's electric charge value, and 
Qc is the minimal energy charge value. Here, each lattice 
site is viewed as a "quantum dot" for which Qc is fixed 
by the (implicit) chemical potential. The last term in the 
Hamiltonian of Eq. (6) specifies that U is the energy cost 
for promoting two electrons into a vacancy and a dou- 
blon. When U/t is large, one can project out doublons 
and consider a restricted subspace of electrons and va- 
cancies only. The local U interaction can then be taken 
into account in second-order perturbation, as shown in 
Fig. 6, renormalizing the coupling constant of the spin 
exchange term to J = Jq + 4t^//7, i.e. the magnon me- 
diated interaction. 



C. Hilbert space truncation for anyons 

For anyons we proceed in the same way as for the elec- 
tronic Hubbard model to derive a simpler effective model 
of the low lying states, assuming that the charging energy 
is the largest energy scale and can be integrated out. In 
Figs. 5(b,c), we model the energy costs for a given quasi- 
particle type using a quantum dot, which again has a 
quadratic Coulombic charging energy Fcoui(^), but also 
has an energy shift 



(8) 



where is the conformal scaling dimension correspond- 
ing to topological charge a, which depends on the (neu- 
tral) topological charge of the quasiparticle^^. We note 
that conformal scaling dimension A^ is the sum of the 
left and right conformal weights A^ = ha 

For large charging energy we can restrict ourselves to a 
low-energy subspace that contains only two quasiparticle 
types, as indicated in Figs. 5(b,c). In the case of the Ising 
anyon chain, we only allow the lowest energy quasiparti- 
cles (/, 0) and (cr, e/4) to be localized on a given site. The 
quasiparticles (/, e/2) and correspond physically 
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to double occupancies of the quasiparticle (cr, e/4), and 
thus involve a prohibitively large Hubbard- like charg- 
ing energy. The neutral fermion quasiparticle (?/^, 0) is 
gapped, because of the energy associated with the 
mode, but it will be present in our model in the form 
of virtual tunneling processes. To make the quasipar- 
ticles (/, 0) and (a, e/4) nearly degenerate, one has to 
introduce a chemical potential, which combined with the 
charging energy gives the sought after near-degeneracy. 
Similarly, in the case of Fibonacci anyons, we can also 
allow only the lowest energy quasiparticles (V^o^O) and 
(r2, e/5) to be localized on a given site. In this case, the 
neutral Fibonacci quasiparticle (ro,0) is gapped, but will 
be present in our model in the form of virtual tunneling 
processes. In Appendix B, we explain in more detail how 
one can combine the effects of a gate and of the charging 
energy in order to obtain a low-energy sector containing 
two degenerate states, separated from the other excita- 
tions by a gap. 

We note that the fundamental quasiparticles of the MR 
or RR quantum Hall states are described as a product 
of an Ising or Fib anyon model with an Abelian anyon 
model (see Appendix A) which can be associated to the 
electric charge and is, therefore, "additive." Hence, the 
electric charges of the quasiparticles of the relevant sub- 
space need not to be specified any more, because each 
quasiparticle has the same charge. Also, we can drop 
the subscripts for the two lowest-energy Fibonacci quasi- 
particles, i.e. we identify iIjq ^ I and T2 ^ t since only 
one species of r-anyons is involved in the low-energy sub- 
space. In other words, we end up with only one type of 
cr- anyon or r- anyon allowed on the sites and the sites left 
empty are filled with trivial quasiparticles / or vacancies. 
The details of the above mentioned identification are not 
important here, but will be given in Appendix B. 

A pictorial representation of such chains was shown in 
Fig. 1. The charge degrees of freedom can therefore be 
thought of as living on the sites (the localized nontrivial 
quasiparticles carry elementary e/4 or e/5 electric charge, 
in contrast to the vacancies) while the bond variables Xi 
are encoding the anyonic (or spin) degrees of freedom. 
We recall that ax / = a for any quasiparticle a = I^r^a^tlj 
so that the anyonic "spin" is conserved along the empty 
segments of the chain (i.e. with vacancies on the sites). 



D. Interaction between nearest-neighbor anyons 

Let us consider putting N anyons of type <j or r on an 
L-site chain with periodic boundary conditions (PBC), 
i.e. on a closed ring (situated on a torus). When two 
charged anyons sit on nearest-neighbor (NN) sites they 
experience an ordinary Coulomb repulsion V. In addi- 
tion, they interact via an effective exchange interaction 
of magnitude J, which can be derived as in the electronic 
Hubbard chain. For this, we use the (unitary) F-symbol 
transformation shown in Fig. 7, which is a change of ba- 
sis between different fusion tree representations of the 



a 




o 

X 



b 





a 



X Z 



FIG. 7: Change of basis involving the fusion channel of two 
neighboring electrons or anyons and the F-symbol. 
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(b) 
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FIG. 8: Matrix elements describing (exchange) interactions 
between nearest neighbors of (a) electrons, (b) Ising anyons, 
and (c) Fibonacci anyons. In (b), when x — z — y and 
y' can take the values I or ip, making h^'^'^ a 2 x 2 matrix, 
while when neither x nor z equals cr, then y — y' — making 
/ij''^''^ a 1 X 1 matrix. 



states. When we apply the F-symbol of the NN anyons, 
which have charges a and 6, it provides a change of basis 
from the fusion chain basis of Fig. 1 (which we use to 
encode states) to one in which the fusion channel of this 
NN pair of anyons is manifest, as indicated in Fig. 7 by 
the charge label w. 

By analogy with the electronic Heisenberg interaction, 
the exchange interaction between two neighboring anyons 
is given by —JV^^ which favors the vacuum fusion chan- 
nel / for these two anyons. The action of the correspond- 
ing exchange processes on the local fusion tree basis ele- 
ments are shown schematically in Figs. 8(a-c). By using 
the F-symbol change of basis of Fig. 7, the local Hamil- 
tonian elements /i^^^J'^^''^''^^ can be derived, depending 
on (and labelled by) the variables and x^+i on the 
two outer bonds connected to the two NN sites and act- 
ing upon the local spin xi of the inner bond, as shown 
in Fig. 8(a-c). The label denotes the type of anyon 
localized at site i. Severe local constraints greatly reduce 
the number of possible non-zero matrices and matrix el- 
ements, which we give explicitly below. 

Let us first start with the case of two NN spin- 1/2 (lo- 
calized) electrons experiencing an AF exchange interac- 
tion, i.e. for which the fusion outcome in the singlet chan- 
nel is favored over the triplet channel. In the usual spin- 
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basis, this is just the Heisenberg term — J(l/4 — -S^+i). 
However, we work in the fusion chain basis, as pictured 
in Fig. 1(a). Thus, we need to know the F-symbols de- 
scribing the change of basis as given in Fig. 7, for the 
case of SU{2) spin-l/2's, i.e. = a^+i = 1/2. The 
F-symbols are closely related to the Wigner 6j-symbols 
[see, for instance Ref. 49 or the 811(2)/^ F-symbols with 
q = 1 (i.e. k = oo) in Appendix A]. The F-symbols of 
interest here are given by 



■^i — 1 5 2 ' 2 



(-1) 



Xi-i-\-Xi + i-\-l 



i =0, 

lis j i23 J 



(9) 



where = 0, 1 is the total spin of the two spin-^'s, and 
denote the 6j-symbols. 



In particular, if x^+i = Xi 
and Xi are fixed to be Xi = Xj 



-1 ± 1, then the value of Xi 
_i ± ^ and Xi = 1, and the 
resulting F-symbol is just a number, namely 



-1±1 



= 1. 



(10) 



In the case that = x^+i = 0, we must have Xi = 
1/2 and Xi = 0. The associated F-symbol is again 1, 



1. 



(11) 



The only case for which the F-symbol has rank two 



is when Xj 



_i = Xi^i = s > ^, givmg Xi 
Xi = 0,1. The F-symbol takes the from 



5 ± ^ and 



F 



1 1 



1 



(12) 



0, and the 



where the first column corresponds to x 
second one x = 1. 

With the knowledge of the F-symbols, we can con- 
struct the Hamiltonian [see Fig. 8(a)], which symbolically 
takes the form 



iXi- 



J 



f: 



Xi,0 



(13) 



where we included the coulomb interaction V, because 
the electrons occupy neighboring sites, and we favor the 
spin-0 channel (implicitly it is only non-zero if the dia- 
gram is allowed by the fusion rules). Explicitly, we find 

that hxi 



-1,1/2,1/2, 



Xi—\ ^i+l 

s > 0, we have 



V in the case that x^+i 7^ For 
V — J , and for 



0, we have /iq'"*^^^'"*^^^ 



/,.,l/2,l/2 



V 

2s+l 



Js Jy/aCl+a) 
28+1 2s+l 

J(5 + l) 

2s+l 



V 



(14) 



The Hamiltonian in the case of Ising and Fibonacci 
anyons [see Figs. 8(b,c)] is obtained in the same way as 
we did above for spin-1/2 electrons. The most important 
necessary ingredient are the F-symbols, which can be 
found in Appendix A. 

For Ising anyons, the non-zero matrix elements are lim- 
ited to 



I ,(J,a 














a, a 


I ,(7, a 








a, a 





= V-J 

= v 



and 



V-J/2 -J/2 ■ 
-J/2 V - J/2 



(15) 
(16) 

(17) 



where the basis used to write the matrix is {/, ?/^}. 

The non-zero matrix elements of the Fibonacci chain 
are given by 



/,r,r 



V-J 



and 



'V-J/<P -J/(/.3/2- 
-J/<^3/2 V-J/<j)^ 



(18) 
(19) 

(20) 



where (j) is the golden ratio and the matrix is written in 
the basis {/, r}. 



E. Anyon "hopping" 

Finally, we have to consider the possibility for quasi- 
particles (including the vacuum /) to move on the lattice 
and gain kinetic energy. The (effective) physical hop- 
ping processes are shown in Fig. 9. In such a move, 
the entire quasiparticle, including the electric charge and 
spin/topological charge, is transported from one lattice 
site to a vacant site that is adjacent to it. (Generally, 
"hopping" may involve transfer of a quasiparticle to an 
occupied site, but we do not consider such processes in 
our models.) Note that the magnitude of the hopping t 
is not affected by the truncation of the Hilbert space to 
the reduced space of the low-energy quasiparticles. Note 
also that the sign of t is irrelevant, so one can assume 
t > for simplicity. 

When I J| is large in comparison to t and V, the system 
phase-separates i.e. the anyons tend to form large clus- 
ters of higher density p ~ 1. In contrast, for larger kinetic 
energy (i.e. t) and/or repulsion between the anyons, the 
system remains homogeneous. This is the regime of in- 
terest here and, for p = 2/3 and p = 1/2, we have found 
that it is realized for t = | J|, even when V = 0, or larger 
t/\J\ values. Note that for convenience we assume V = 
throughout and we have explored the models for values 
of J ranging from —t to t. 
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Hardcore bosons 



(a) 



(b) 



(c) 



Xi X2 X2 



,X'l ,X'i ,T2 



FIG. 9: Tunneling process (or "hopping") of (a) an electron, 
(b) an Ising anyon, and (c) a Fibonacci anyon. 



IV. CHARGE DEGREES OF FREEDOM FOR 

J = 

This section is a "warm-up" for the real itinerant anyon 
models, starting with the simple example of identical 
bosons, and describes how making them distinguishable 
introduces a twist in the periodic boundary conditions. It 
therefore explains the J = part of the spectra, without 
the non-Abelian complications. 



A. Hard-core bosons 

We start with a periodic chain of size L filled with N 
bosons. We consider the case where two bosons cannot 
occupy the same site due e.g. to an infinite on-site repul- 
sion (hard-core constraint). Such a system of hard-core 
bosons (HCBs) can be mapped via a Holstein-Primakoff 
and Jordan- Wigner transformation onto a gas of spin- 
less fermions. In ID, the effect on the spectrum due 
to the difference in statistics can simply be accounted 
for by adding to the fermions an extra phase shift of tt 
through the ring (when the particle number N = pL is 
even). Therefore, the HCB many-body spectrum can be 
obtained by filling up N states of a (fermionic) cosine 
band, 



^hcb(p) 



-2t cos 

jip) 



(21) 



where {j{p)} is a set (labelled by p) of an even number N 
of different integers and the momenta are all shifted by 
7t/L. The spectrum (for t = 1) is displayed in Fig. 10(a). 
As expected, the HCB spectra exhibit linear quasiparticle 
dispersions centered at momenta K = and K = Kc = 
27rp [or 27r(l - p) for p > 1/2]. 

For later use in the case of anyons, it is of interest to 
introduce an external magnetic fiux ^ext or equivalently 




momentum K 



K+P<^e> 



FIG. 10: Spectra of 16 (hardcore) bosons moving on a 24- 
site chain with PEG. (a) Spectrum for zero magnetic flux 
through the ring. The linear dispersions vs momentum K are 
shown at K = and K = zbiCc. (b) Spectrum vs. external 
(continuous) magnetic flux (/)ext (multiplied by the density p). 
Each discrete level (black dots) leads to a (parabolic) branch 
of excitations. 



an Abelian U(l) flux through the ring. The new spec- 
trum ^hcb(p, 0ext) depends now on both discrete and 
continuous variables p and (/)ext ? 



EuCb{P, ^ext) 



-2t ^ cos 

jip) 



27r . 1 . 0ext 



, (22) 



plotted in Fig. 10(b) (for t = l).The state labelled by 
(p, ^ext) now carries an arbitrary (continuous) total mo- 



mentum K = Kn 



with Kn 



27r 
L 



B. Fermions and anyons 

We now make the (quasi)particles distinguishable^ i.e. 
we introduce some internal degrees of freedom which can 
be e.g. the spin-1/2 components of the electrons or the 
anyonic degrees of freedom of the Ising a or Fibonacci r 
anyons. The resulting model is the same as considering 
the above t-J models in the limit where J = (t can be 
set to 1). This limit, where the energy scale of the anyonic 
degrees of freedom are set to zero, is of great interest since 
it provides insight on the nature of the charge excitations. 

In that limit, the mobile anyons are expected to still 
behave as HCB. However, the extra internal degrees of 
freedom (with zero energy scale) should provide extra 
features on top of the HCB spectrum. Because of the 
anyonic (or spin) degrees of freedom linked to them the 
charged bosons are no longer indistinguishable particles 
and, on a torus or a ring, hopping of a particle across the 
"boundary" cyclically translates the labels of the fusion 
tree. To recover the same labeling, in general, all N 
particles must be translated over the boundary. Thus, 
one distinguishable particle hopping over the boundary 
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K+pc|)g^^ momentum K K+p(j) 



FIG. 11: Low-energy spectrum of a 24-site t-J chain at J = 
and density p = 2/3, obtained numerically using Ising anyons. 
Yellow circles correspond to the spectrum at (j)e^t = as a 
function of momentum K. The "parent" charge excitations 
(HCB at (/)ext = 0) are shown by full (black) circles. Adding 
an external flux (/)ext through the ring is equivalent to shift 
K by p(/)ext : the (red) x symbols correspond to the spec- 
trum at (/)ext = 7r/4 restricting K G [— 7r/3,7r/6] (so that 
K e [-7r/6,7r/3]). 



FIG. 12: Low-energy spectra vs. momentum K oi a, 18-site 
t-J chain with J = at anyon densities (a) p = 1/2 and 
(b) p = 2/3. The notations here are the same as in Fig. 11, 
but these results are obtained numerically using Fibonacci 
anyons. The low-energy spectra of the 18-site HCB chain at 
the same densities (parent excitations) are shown by black 
dots. Data for an external flux (/)ext = 7r/5 are shown in (b). 
The minimum of the spectrum occurs at momentum K = 
or K = TV, depending on the parity of the number N = pL oi 
quasiparticles. 



has the same effect as a phase shift (pn = ^tt-^ (with n an 
integer). Hence the complete J = electronic/anyonic 
spectrum (at zero external flux) is given by the union of 
all HCB spectra taken at all discrete values of 0^, 



charge 



(23) 



A momentum shift is induced by the U(l) flux, given by 
p(j)n i.e. 27r^, an integer multiple of The states then 
carry (discrete) total momenta, 



(24) 



where Kp are the momenta of the HCB eigenstates at 
^ext = 0- For convenience one can distribute the phase 
shift equally on the bonds to preserve translational in- 
var iance and one gets. 



"^charge 



-2t ^ cos 



27r 



(25) 



Spectra for Ising and Fibonacci anyons obtained by ex- 
act diagonalization (ED) for J = are shown in Fig. 11 
and Fig. 12 (a, b). As expected, Eq. (25) matches ex- 
actly the numerical exact diagonalization results. From 
the above considerations, it is then clear that the eigen- 
energies lie exactly on top of the parabolas corresponding 
to the "optical" excitations of the HCB. In other words, 
each state in the HCB spectrum is extended into a dis- 
crete set of levels on a parabola - the same parabola that 
one gets by adding flux (an Abelian phase), as checked 
numerically. 



It is important to notice that the J = spectrum does 
not depend on the internal degrees of freedom and, hence, 
on the nature of the quasiparticles, i.e. whether they are 
electrons, Ising anyons, Fibonacci anyons, or distinguish- 
able bosons. However, the respective Hilbert spaces are 
very different which means that the corresponding eigen- 
functions and degeneracies differ completely. In addition, 
the way the very large degeneracy of each level is lifted 
by any finite exchange interaction (see Fig. 14 discussed 
later) depends crucially on the type of particles. 



C. External magnetic flux 

When the anyons experience an arbitrary external 
flux 0ext, the above formula can be generalized to 
^chrrge(^ext) = £^hcb(p, (/>n + (/>ext). It then bccomcs ap- 
parent that the J = energy spectrum does not depend 
on the momentum Kp^n external flux ^ext separately 
but rather only on the "pseudo-momentum" combination 



/)ext- Hence, one can define a spectrum de- 



pending on both discrete and continuous variables, 

£^charge(p, ^) = ^HCb(P, ^) • 



(26) 



The curvature of the ground-state energy 
S^£^charge(0, <^>)/9^^ is directly proportional to the 
optical (Drude) weight quantifying the potential of this 
system to conduct. 
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FIG. 13: Energy spectrum of a dense Ising-anyon chain (p = 
1) of length La = 16. The ground state energy has been 
subtracted. 



V. DENSE ANYON MODELS 

To complete the warmup to describe the full anyonic 
t-J model, we briefly discuss the dense anyon models at 
p = 1. These models have precisely one anyon per site, 
which are, hence, immobile due to the hard-core con- 
straint. Every anyon interacts with its two neighbors, 
and only the sign of the interaction strength J is rele- 
vant. These models, introduced in Ref. 50, are the any- 
onic versions of the Heisenberg spin chains. We will only 
consider the spin- 1/2 versions in this paper. 

In the case where only a nearest-neighbor two-body in- 
teraction is present, the spin- 1/2 anyonic chains are all 
critical, and their energy spectra are described by well 
known conformal field theories. Starting with the Ising 
anyons, we note that, due to the fusion rules, the degrees 
of freedom on the fusion chain are forced to form a pat- 
tern of alternating frozen a bonds and bonds fluctuating 
between / and ip. For these later bond variables, the in- 
teractions of Fig. 7(c) are exactly those of a critical Ising 
model in transverse field whose corresponding GET has 
central charge c = 1/2. This is irrespective of the overall 
sign of the interaction, although the momenta at which 
the various states occur differ depending on the sign of 
J. 

In the case of Fibonacci anyons, changing the sign of 
J does alter the critical behavior of the chain. In the 
case of anti-ferromagnetic interactions (favoring the triv- 
ial fusion channel of two neighboring Fibonacci anyons), 
the critical behavior is described in terms of the c = 7/10 
tri-critical Ising model, with low-lying, linearly dispers- 
ing modes occurring at momenta K = and K = ir. 
For ferromagnetic interactions, the critical behavior is 
instead described by the c = 4/5 3-state Potts model, 
which exhibits low-lying modes at = 0, 27r/3, 47r/3. 

The behavior described above can be obtained by map- 
ping the models onto exactly solvable two-dimensional 



height models, introduced by Andrews, Baxter and 
Forrester. The result for spin-1/2 anyons associated 
with SU(2)/e is that in case of anti- ferromagnetic inter- 
actions, the critical behavior of the chain is given in 
terms of the /c-critical Ising model, while ferromagnetic 
interactions give rise to the critical behavior of the Zk- 
parafermions. 

We would like to stress that, although one can analyt- 
ically obtain the critical behavior of the dense anyonic 
chains, it is not possible to obtain the energy spectra 
in full detail for finite size systems for /c > 3. To obtain 
these, one must employ numerical techniques, such as ex- 
act diagonalization (see e.g. Fig. 13, for the k = 2 Ising 
case, which can also be obtained exactly). In describing 
the full spectra of the anyonic t-J models, we will make 
use of the spectra of the dense anyonic models described 
here, as obtained from exact diagonalization. 

We will denote the length of the dense anyon chains 
by La- The energies of the dense anyon chains are de- 
noted by £^anyon(^), whcrc the integer m labels the eigen- 
states, which have momenta km that are integer mul- 
tiples of 27T/La. Next, considering chains of length L 
and at anyon densities p < 1, we shall define the cor- 
responding "squeezed chains" of dense anyons of length 
La = pL = A/", in which the vacancies (or trivial quasi- 
particles) have been removed. 



VI. MANY-BODY SPECTRA OF THE 
ANYONIC t-J MODELS 

Having described the spectra of the HCB system, in 
the presence of external flux, as well as the spectra of 
dense anyon models, we are now ready to describe the 
spectra of the full, itinerant anyon models. We will label 
the various energies as Ep^rn-, where the labels p and m 
refer to the (renormalized) HCB spectrum and the dense 
anyon chain, respectively. We will also explain the subtle 
coupling of the momenta. 



Separation of charge and anyonic degrees of 
freedom 



We now consider to the full J > spectra of the itin- 
erant models. We solve the effective anyonic models on 
small periodic rings using exact diaginalization. For Ising 
anyons, a 24-site chain is studied at density p = 2/3 
{N = 16 quasiparticles). For Fibonacci anyons (which 
have a larger Hilbert space) we consider a 18-site chain 
at density p=l/2(A^ = 9 quasiparticles) and p = 2/3 
(TV = 12 quasiparticles). We choose \J\ < t for which 
the system remains uniform and does not phase separate 
(which occurs for larger J). Note that the sign of J is ir- 
relevant for the energy levels in the case of Ising anyons, 
though the momenta at which the various states occur 
differs depending on the sign of J. 
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FIG. 14: A zoom-in on the low-energy spectra vs. K oi a, 24- 
site Ising anyon t-J chains at density p = 2/3 for (a) J = and 
(b) a small value of J/t = tan(7r/50) ^ 0.06. The Lanczos 
algorithm with 800 iterations is used so that, in the shown 
energy window, most of the eigen-energies have converged to 
within a relative error of 10"^^ (a few not-fully convergea 
levels are not shown). The (red) -h symbols correspond to 
the sum of the (computed) lowest charge branch (continuous 
red line) with all the expected anyonic excitations. (See text 
for more details.) 



FIG. 15: Low-energy spectra vs. K oi a, 24-site Ising t-J 
chain at density p = 2/3 for (a) J = and (b) J = t = 1/a/2. 
The data at (/)ext = are shown by yellow circles. The parent 
charge excitations are shown by black dots ((/)ext = 0) and, 
as a function of pseudo- momentum K = K -\- (/)extp (varying 
0ext), by continuous lines of different colors. The + and x 
symbols added for comparison correspond to the sum of the 
charge and expected anyonic excitation spectra. (See text for 
more details.) The colors of these symbols are the same as 
their parent charge excitation parabolas. 



The low-energy spectra (|J| = t = 1/a/2) of the itin- 
erant Ising and Fibonacci anyonic chains are shown in 
Figs. 15(b) and Fig. 16(a,b), respectively. These seem 
very different from the J = limit studied above (and 
shown again in Fig. 15(a) for comparison). To under- 
stand such spectra, let us first consider a zoom-in on the 
low-energy region and compare the spectra at J = and 
at a small value of J, as shown in Figs. 14(a,b). This re- 
veals that the highly degenerate J = charge excitation 
parabola is being split by the magnetic interaction into a 
complex spectrum with a spread in energy proportional 
to JL. When J ~ t/L^, the spectra originating from each 
parabola start to overlap as expected in Figs. 15(b) and 
16(a,b). Despite the apparent complexity of the J ^ 
spectrum, we shall be able to express all excitations as 
the sum of an anyonic excitation and a charge excitation, 
extending the concept of spin-charge separation familiar 
for ID correlated electrons to the case of a ID anyonic 
interacting system. To complete this task, we first estab- 
lish from simple considerations the "recipes" to construct 
separately the expected charge and anyonic spectra. In 
a second step, we show how the numerical spectra of the 
t-J anyonic chains can be seen as a subtle combination 
of the above two spectra. 

The Bethe ansatz results^^^^ for the J ^ elec- 
tronic t-J chain suggest that the anyonic contributions, 
^anyon(^), to the cxcitatiou spectrum of the itinerant 
anyon chain are those of the squeezed periodic chain of 
localized anyons produced by removing all vacant sites, 
which has the resulting length La = N = pL. Here, 



the integer m labels the eigenstates of momenta km, 
which are multiples of 27r/ La. Such a spectrum can be 
computed separately by exact diagonalization and agrees 
very well with the CFT predictions, even on small chains 
{La = 12, 16). In particular, it shows a (linear) zero en- 
ergy mode at zero momentum and at a characteristic mo- 
mentum ka, where ka = tt for Ising and J > Fibonacci 
chains and /c^ = 27r/3 for J < Fibonacci chains. The 
coupling constant providing the scale of the anyon spec- 
trum is expected to be weakly renormalized from J to 
J J in the doped system, where 7 is a factor of order 1 
that is to be adjusted as we described below. 

To construct the expected charge excitation spectrum 
at finite J, we use our understanding of the charge ex- 
citations in the J = limit. Starting from J = 
and turning on J gradually, one can a priori adiabati- 
cally follow the original (<l> = 0) HCB excitations evolv- 
ing in, what we call, the parent charge excitations at 
J 7^ (labelled by the same integers p and at the same 
momenta Kp). As for J = 0, changing the momen- 
tum K of di charge excitation amounts to introducing 
a total phase shift (or flux) ^ = K/p. Hence, by in- 
troducing "twisted boundary conditions" one can com- 
pute numerically the (almost parabolic) branch of exci- 
tations ^charge (^7 ^) associatcd to each parent excitation 
(labelled by p). Note that each branch is "renormal- 
ized" by J so that, strictly speaking, the charge spec- 
trum is no longer associated to non-interacting spinless 
fermions (i.e. to HCBs) and, hence, is no longer given 
by a simple analytic expression as in Eq. (26). How- 
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Fibonacci anyons 




FIG. 16: Low-energy spectra vs. K oi a, 18-site Fibonacci t-J 
chain at density p = 2/3 and \J\ = t = 1/a/2 for both (a) J < 
and (b) J > 0. The data at (/)ext = are shown by yellow 
circles. The parent charge excitations are shown by black dots 
(0ext = 0) and, as a function of pseudo-momentum K = K -\- 
P0ext (varying 0ext), by continuous lines of different colors. 
The + and x symbols added for comparison correspond to the 
sum of the charge and (expected) anyonic excitation spectra 
(see text). The color of these symbols is the same as their 
parent charge excitation parabola. 



ever, if different parent states lie on the same branch, 
they are still exactly spaced apart by integer multiples of 
= 27r i.e. ^charge (p, ^) = ^charge (p', ^ + 2fc7r) whcrc 
Kp' — Kp = 2/cp7r, for some integer k. 

We now explain how to construct the full excitation 
spectrum by simply considering that (i) the charge de- 
grees of freedom are subject to a phase shift in the bound- 
ary conditions and (ii) the anyonic degrees of freedom are 
the ones of the squeezed periodic anyonic chain. Accord- 
ing to the above arguments, the energy excitation spec- 
trum should be given by adding the two contributions, 

-^p,m — -E^chargei^Pi ^m) H~ -^anyon(^) • (^'^j 

A natural prescription is to simply add the momenta: 
K = Kp-\-pkm' In other words, we assume that the phase 
shift experienced by the charged "holons" coincides with 
the total momentum km = (where rim is an inte- 

ger) of the anyonic eigenstates defined on the squeezed 
(undoped) chain. These rules for adding charge and any- 
onic momenta are therefore assumed to be similar to the 
J ^ Bethe ansatz. 

We now wish to verify that proper assignments of 
the true energy levels according to the form given by 
Eq. (27) can indeed be made accurately. First, we con- 
sider adding a very small exchange coupling J which lifts 
the very large degeneracy of the low-energy parabola of 
the HCBs (see Fig. 14). For finite J the first charge 
branch (with £^anyon = 0) originating from the zero- 
momentum ground-state of the model (which we assign 
p = 0) can be computed by adding an Abelian flux to 



FIG. 17: Renormalization parameter 7 of the energy scale 
of the anyonic degrees of freedom for Ising and Fibonacci at 
p = 2/3, computed on L = 24 and L = 18 chains, respectively. 



the system. It is then possible to construct the expected 
set of combined charge plus anyon excitations £^o,m by 
adjusting the renormalization factor 7 to get the best fit 
to the exact low-energy levels. Although there is only 
one free parameter, it is remarkable that all anyon exci- 
tations above the lowest charge parabola can be assigned 
very accurately as seen in Fig. 14(b). 

When I J| ~ t, charge and anyonic excitations have the 
same energy scale and one must proceed step by step, 
sequentially constructing the sets of levels correspond- 
ing to increasing charge index p. The two "secondary" 
parent charge excitations corresponding to exact eigen- 
states of the system with momenta Kp = {p = 1,2) 
lie on the same p = branch, as seen on Figs. 15(b) 
and 16 (a, b). These states lead to the secondary level of 
combined excitations Ei^m and £^2,m7 with no further ad- 
justable parameter. Recall that Kc = 27rp if p < 1/2 and 
Kc = 27r(l — p) for p > 1/2. Next, in a second step, 
we identify the lowest not yet assigned excitations at mo- 
menta Kp = ±27r/L as the subsequent pure charge exci- 
tations (and assign them the labels p = 3,4). Following 
these levels adiabatically under the addition of a flux en- 
ables us to construct the corresponding charge branches 
and locate the secondary pure charge excitation at mo- 
menta Kp = ±(27r/L + Kc) (called p = 5,6). Then, as 
before, the combined excitations Ep^mi P = 3,..., 6, can 
be constructed. One can repeat this procedure (going 
up in energy) until the level density and the number of 
charge branch crossings becomes too large to make pre- 
cise assignments. In practice, we have identified up to 
p = 11 pure charge excitations and their corresponding 
low-energy combined anyonic- charge excitations for the 
Ising chain, as shown in Fig. 15(b). For the Fibonacci 
chains, we have identified up to p = 9 pure charge ex- 
citations and their corresponding low-energy combined 
anyonic-charge excitations, as shown in Fig. 16(a,b). 

Our results show that the anyonic energy spectrum 
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is basically given by the same type of Bethe ansatz as 
for the electronic t-J model (in the J /t limit)^^'^^ 
and, in particular, (i) the J = charge excitation spec- 
trum is exactly the same^ (ii) the spin excitations also 
correspond to the squeezed localized chain, and (iii) the 
rules for adding charge (holon) and spin/anyon momenta 
are identical. In addition, the numerical spectra agree 
very well with the sum of the spin (provided some renor- 
malization 7 of the energy scale, as shown in Fig. 17) 
and charge spectra (constructed independently) with the 
above-mentioned rule for momentum conservation. We 
believe the small deviations can be attributed to finite 
size effects (which vanish when J/t 0). Interestingly, 
the J ^ limit of the renormalization parameter 7, 
7(0) 0.5515 for /) = 2/3, is independent on the anyon 
type (as the J = charge excitation spectrum). In fact, 
in this limit, it should only depend on the probability of 
having two neighboring anyons. 



hold in the case of all low-energy eigenstates of ID non- 
Abelian anyons at J 7^ 0, up to finite size corrections. We 
speculate that the eigenfunctions can be approximately 
written as, 

Vn 

where yj are the position of the (site) anyons on the L-site 
chain and Xi are the bond variables associated to them 
(see Fig. 1). Here, ^charge eigenstates (labelled by 

p) of an interacting L-site HCB chain in the presence of a 
twist (i.e. fiux) km in the boundary conditions and X^yon 
are the eigenstates (labelled by m) with momentum km 
of the interacting (undoped) anyonic chain of La = pL 
sites. 



VII. DMRG STUDY OF THE MODEL 



B. Anyonic and charge collective excitations 

Due to the above decoupling, the collective any- 
onic and charge excitations can be deduced easily from 
the above excitation spectra, by just applying selection 
rules. Charge excitations occur between different charge 
branches at constant m with energy transfer 



-E'charge(P ? ^m) -^charge (P 7 ^m) 7 (^^) 



and momentum transfer K = Kpi — Kp. Anyonic exci- 
tations are characterized by Ap = and are then given 

by 

^p]m' ,m -E'anyon('^ ) -^anyon(^) 

H~ -E^charge(P7 ^m') -^chargei^Pi ^m) i (^9) 

with momentum transfer K = p{km' — km)- Note that 
the last two terms in Eq. (29) give a finite size correc- 
tion in the energy of order 1/L^. In the thermodynamic 
limit, zero-energy anyonic excitations occur at momen- 
tum Ka = pka (and 2Ka if different), where ka is the 
characteristic momentum (introduced above) of the zero- 
energy mode of the pure chain. The location of both 
charge and anyonic zero-energy modes are indicated in 
Fig. 15(b) for Ising anyons and in Figs. 16 (a, b) for Fi- 
bonacci anyons. 



C. Form of the eigenstates 

We now discuss briefiy the structure of the eigenfunc- 
tions. Eq. (2.14) of Ref. 52 established that the ground- 
state of the J ^ limit of the electronic t-J chain can be 
written exactly as a product of a charge HCB wavefunc- 
tion times a spin wavefunction identical to the ground- 
state of the ID = 1/2 Heisenberg model. Our results 
suggest that a similar product structure might in fact also 



We can use DMRG to compute the resulting CFT cen- 
tral charge from the analysis of the von Neumann entan- 
glement entropy (EE) of an open chain (of length L) cut 
into two subsystems'^. We focus here on the case of the 
diluted Fibonacci t-J chain with p = 2/3 and J > 0, for 
which we expect the anyonic part to be described by a 
c = 7/10 tri-critical Ising CFT. In Fig. 18, we plot the 
ground state's EE 



Sa = Sb = -Tr [pA ^og Pa] 



(31) 



between subsystems A and which are two connected 
segments of the open chain, as a function of the position 
of the cut along the chain. The calculation is somewhat 
non-standard (compared to usual spin systems) because 
the anyonic fusion tree bond variable Xj labeling the jth 
link, which connects the two subsystems, i.e. the link 
across which one "cuts" the system in two, is shared by 
both subsystems of the chain. This shared link variable 
characterizes the overall topological charge of each sub- 
system. The reduced density matrices 



pA Tr^p, pB Tr^p 



(32) 



of subsystems A and B are block diagonal with respect 
to this variable, i.e. 



PA = ^PxjpA,Xj, 



(33) 



where Px^ = Tr[IIa;^.p] is the probability that the state will 
have topological charge xj on the jth link, and pA,xj = 
^TtaI^XjP] is the reduced density matrix for subsystem 
A after projecting the jth link's variable onto the value 

Xj. 

In order to verify that this Fibonacci t-J chain results 
in a CFT with central charge c = 1.7, we fit to the for- 
mula^'^^ 



= a^bpxj=r + - log 



L sin 



7TJ 



(34) 



15 



C/5 
>^ 
Q. 
O 



1.4 




■ ■ . . ^DMRG - 






ii w W w VAa a f '^'^ 


1.2 






1.0 






0.8 








& \ 


0.6 




Q _ 

1 1 1 1 1 1 1 



10 20 30 40 50 60 70 

position X 



FIG. 18: Entanglement entropy obtained using DMRG for an 
open chain of length L = 72, with N = 48 Fibonacci anyons, 
and J/t = 0.3. The EE predicted for a CFT with central 
charge c = 1.7 is plotted and the agreement is seen to be 
excellent. 



where a and b are fitting parameters. The first term is a 
non- universal constant, which can include universal con- 
tributions, such as a boundary entropy. The second term 
is a phenomenologically motivated correction that is pro- 
portional to a local kinetic energy, i.e. Pxj=T = (^iink(j)) 
where niink(j) is the density (occupation) operator of the 
jth link, and can include a contribution due to the bound- 
ary between the two subsystems. The third term is de- 
rived from CFT. We find the best fit for the parameter 
values a = 0.31185 and b = -0.35547. As seen in Fig. 18, 
the agreement between the numerical results and the val- 
ues provided by this expression is excellent. 



VIII. CONCLUSION AND OUTLOOK 

Motivated by the possible realization of non-Abelian 
Ising and Fibonacci quasiparticles in quantum Hall states 
and Majorana heterostructures and the importance of 
understanding their edge modes, we have investigated 
what happens if itinerant and interacting (charged) non- 
Abelian anyons are confined on a ID chain, subject to a 
strong charging energy. Following a standard procedure 
for strongly correlated electronic systems, we have con- 
structed simple low-energy effective models by truncat- 
ing the Hilbert space to the relevant low-energy particles. 
Integrating out the high-energy virtual processes yields 
an "exchange" interaction between anyons, which phys- 
ically favors a particular fusion channel. The effective 
model generically takes the form of an anyonic t-J model, 
containing the exchange interaction (J) and the rate (t) 
of anyon "hopping" between nearest-neighbor sites. The 
central result of our work is that anyons fractionalize into 
their charge and (neutral) anyonic degrees of freedom. 
This phenomenon closely resembles and generalizes the 
well-known spin-charge separation in electronic Luttinger 
liquids. Incidentally, the numerically verification based 
on the identification of the many-body levels turned out 
to be more transparent for anyons, due to the absence 



of marginally irrelevant operators in the field theory de- 
scription. The anyon fractionalization justifies a posteri- 
ori the treatment of the edge theories of these topologi- 
cal phases as a direct product of the charge and neutral 
non-Abelian modes, even though the electric charge is 
not localized in current setups. 

We note that the ID electronic t-J model exhibits an 
exact supersymmetric point^^'^"^ at which the full excita- 
tion spectrum can be obtained using the Bethe ansatz^^. 
It is left for future studies to investigate whether such an 
integrable point also exists in ID anyonic t-J models. 

Our simple description of interacting itinerant anyons 
now enables the investigation of realistic setups for ma- 
nipulating and/or braiding anyons for future quantum 
computation. It is also easy to extend this study to quasi- 
ID systems with more than a single conduction chan- 
nel: anyonic t-J "ladders" could mimic such a case, fol- 
lowing the procedure for localized non-Abelian anyons^^. 
Whether fractionalization survives in two spatial dimen- 
sions is another important issue. Localized anyons were 
shown to nucleate into a novel gapped quantum liquid^^ in 
two dimensions and a similar scenario might take place 
for itinerant anyons, with e.g. the anyonic degrees of 
freedom becoming gapped. 
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Appendix A: Examples of Anyon Models 

In this Appendix, we give detailed descriptions of the 
Ising, Fibonacci, and SU(2)/e anyons models, and explain 
where they occur in non-Abelian quantum Hall states. 

1. Ising anyons 

The Ising anyon model is derived from the CFT that 
describes the Ising model at criticality^^. It is related 
to SU(2)2 as its CFT can be obtained using the coset 
construction SU(2)2/U(1)4. It has topological charges 
C = {/, cr, i/j} (which respectively correspond to vacuum, 
spin, and Majorana fermions in the CFT, and are some- 
times denoted 0, ^, and 1, because of the relation with 
SU(2)2). The anyon model is described by (listing only 
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the non-trivial F-symbols and i?-symbols, i.e. those not 
Usted are equal to one if their vertices are permitted by 
fusion, and equal to zero if they are not permitted): 



C = {/, cr, ?/^} , / X a = a, a X a = I ^ il)^ 
(J X ip = cr^ ip X ip = I 



\ rpcraa] 



1 1 

1 -1 
V2 V2 J e/ 



RY 



TDCTCT 



dj = = 1, da- = a/2, V 



topological charges C = {/, r} (sometimes denoted and 
1, respectively, because of the relation with SU (2)3) and 
is described by (listing only the non-trivial F-symbols 
and i?-symbols): 



C = {/, r}, 1x1 = 1, I X T = T, r X r = I ^ r 


[^rie/ = 




ef 


RY = e-^'^^/^ RY = e'^""/^ 


di = l, dr = (t), P = v/0 + 2 


Oi = 1, er = e'^ 



where (j) = is the Golden ratio. 



where e, / G {/, i/j}. 



3. SU(2), 



2. Fibonacci anyons 

The Fibonacci anyon model (also known as SO (3) 3, 
since it may be obtained from the SU (2)3 anyon model 
by restricting to integer spins j = 0, 1, though SO (3)^ 
is only allowed for /c = mod 4; as a Chern-Simons or 
WZW theory, it may, more properly, be equated with 
(02)1) is known to be universal for TQC^^. It has two 



The SU(2)^ anyon models (for k an integer) are "g- 

deformed" versions of the usual SU(2) for q = e*^+2, 
which, roughly speaking, means integers n are replaced 

by [n]^ = '^i/2lg-i/2 • These describe SU(2)^ Chern- 
Simons theories^^ and WZW CFTs^^'^^, and give rise 
to the Jones polynomials of knot theory^^. Their braid- 
ing statistics are known to be universal for TQC^^ all 
except A: = 1, 2, and 4. They are described by: 



min{ji+j2,fc-ji-i2} 



c = {o,^,...,|}, jixh= E 



j=\jl-j2\ 



j 



Ji 32 J12 

js j j23 



V ^ ^ [j3 J J23 

A ( ji , j2 , jl2 ) A ( ji2 , j3 , j) A {j2 , js , j23) A ( ji , ^23 , j) 



i-ir[z+i],\ 

-j2-jl2]Az-jl2-j3-j]Az-j2-j3-j23]a-[z-jl-j23-j]a- 



[ jl + J 2 + js + j -Z] + j 1 2 + js +i23 - 2;] ! [j2 +i 1 2 +i + ^23 - 2;] , 



A(ji, j2, js) = 



[- J 1 + J 2 + js ] „ ! [j 1 - J 2 + js ] „ ! [j 1 + J 2 - js 



[jl+j2+js + l]„ 



^ n H, 

m=l 



^ (_l)j-jl-j2 ^I(j(j + l)-jl(jl+l)-j2(j2 + l)) 



^(^) 

sin(^) 



^(fcT2) 



where { is a "g-deformed" version of the usual SU(2) 
6j-symbols (which correspond to (7 = 1), and have been 
calculated in Ref. 70 (see also Ref. 71, for an introduction 
on how to calculate the F-symbols and an implementa- 
tion in Mathematica). The sum in the definition of the 
g'-deformed 6j-symbol is over all integers in the range 



max{ ji + j2 + ji2 ; ji2 + js + j ; j2 + js + j2S ; jl + j2S + j } < 

Z < min{j\+j2+j3+j; j\+j\2+j3+j23; j2+jl2+j+j23}. 
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4. Mo ore- Read, anti-PfafRan, and 
Bonderson-Slingerland Hierarchy States 

The V = 1/m MR states^ are described by a spec- 
trum restriction of the product of the Ising CFT with an 
Abehan U (1). Specificahy, the anyon model is 

MR = Ising xU(l)^lc (Al) 

where the restriction to the anyonic charge spectrum C is 
such that / and ijj Ising charges are paired with integer 
U (1) fluxes, while a Ising charges are paired with half- 
integer U (1) fluxes. The fundamental quasihole of the 
MR state has electric charge e/2m (where the particle 
carries charge — e) and carries Ising topological charge 
(J. The V = 1/2 MR state is a leading candidate for the 
experimentally observed v = b/2 and 7/2 quantum Hall 
plateaus. 

Taking the particle- hole conjugate of the MR state 
yields the aPf state which is another leading candi- 
date for the V = b/2 and 7/2 quantum Hall plateaus. 
The anyon model for the aPf state is simply obtained by 
taking the complex conjugate of the MR state's anyon 
model. 

BS hierarchical states^ may be obtained from the MR 
and aPf states by applying a hierarchical (or, equiva- 
lently, a composite fermion) construction to the U (1) 
sector. The states built on MR may be written as 

BSk= Ising xU(l)^lc (A2) 

where the i^-matrix is determined by the details of the 
hierarchical construction over MR, and the spectrum re- 
striction is similar to before. This produces Ising-type 
candidate states for all other observed second Landau 
level FQH filling fractions (including those observed at 
V = 7/3, 12/5, 8/3, and 14/5). The quasiparticle exci- 
tation spectra of the BS states include excitations that 
carry the a Ising topological charge, but these are gen- 
erally not the unique quasiparticle carrying the minimal 
electric charge. 



5. /c = 3 Read-Rezayi and NASS 

The particle-hole conjugate of the /c = 3, M = 1 
RR state^ is a candidate for v = 12/5, which is con- 
structed from the Za-Parafermion (Pfa) CFT and an 
Abelian U(l). The braiding statistics of this state is 
described by the direct product of anyon models 



RRfc=3,M=i = Pf3 X U (1) = Fib X (A3) 

where the overline indicates complex conjugation and 
Tj^q is an Abelian anyon model (using the notation of 
Ref. 72,73). The fundamental quasiholes of this state 
have electric charge e/5 and Fibonacci topological charge 
r. 



The /c = 2, M = 1 NASS state^o, based on SU(3)^- 
parafermions, is a candidate for v = 4/7. Its braiding 
statistics is described by 

NASSfe^2,M=i = Fib X (Z2) X U (1) X U (1) , (A4) 

where (Z2) is an Abelian theory similar to D (Z2), the 
quantum double of Z2 (a.k.a. the toric code). The two 
U{1) factors describe the charge and spin of the particles. 
Its data is listed in Ref. 72 and also as = 8 in Table 2 of 
Ref. 74. The fundamental quasiholes of this state carry- 
ing Fib topological charge r, and electric charge of either 
e/7 or 2e/7. 

As these theories are the direct product of a Fi- 
bonacci theory with Abelian sectors, the braiding statis- 
tics of quasiparticle excitations carrying the non-trivial 
Fibonacci charge are computationally universal. 



Appendix B: Limiting the number of particles 

In this Appendix, we will provide some details on how 
we can restrict the quasiparticle spectrum to just two 
quasiparticles, namely the trivial "vacuum" quasiparti- 
cle /, and a "fundamental" or "elementary" excitation. 
This excitation has the smallest possible electric charge, 
and is fundamental in the sense that all other excitations 
can be obtained from it by repeated fusion. The reason 
behind this truncation is to come up with a model which 
is tractable, because the Hilbert space grows exponen- 
tially in the number of quasiparticle types. Apart from 
the truncation of the spectrum to two quasiparticle types 
only, we also need that the energy associated with these 
particles is the same. To achieve both goals, one has two 
tools, in principle. In particular, one can consider using a 
gate which couples linearly to the charge, or consider the 
charging energy associated with localizing the anyons on 
quantum dots. As we will show below, by using a gate 
alone, one can arrange the system to have degenerate lev- 
els for the quasiparticles, but it turns out that this lead 
to a degeneracy larger than two. To split this degeneracy, 
the quantum dot charging is essential. 

After we explain how we can restrict the number of 
quasiparticles to these two types, we explain how we can 
map the obtained model to the anyonic models intro- 
duced in Ref. 50 (in the case p = 1). We do not consider 
the MR and /c = 3 RR cases separately, but directly 
consider the general case of the fermionic RR states for 
arbitrary k (which includes MR at /c = 2). In addition, 
we only focus on those aspects which we will need for our 
purposes in this paper. 

We start by decomposing the operators creating the 
different types of particles into two pieces, one associated 
with the non- Abelian statistics, the other with the elec- 
tric charge of the quasiparticles. The operator describ- 
ing the fundamental quasiparticles, which have charge 
e/{k -h 2) is of the from ^le''^/(VHk+2)) ^ ^^^^^ ^1 ^ 
parafermion field, corresponding to a and <Ji for /c = 2 
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and A: = 3, respectively, using the notation of section II B. 

The vertex operator e^'^^^y^^^\ where is a chiral 
bosonic scalar field, gives the charge of the quasiparticle. 

The energy of the quasiparticles in a finite geometry 
(such as the quantum dots used to localize the quasipar- 
ticles) are proportional to the scaling dimensions of the 
fields creating the particles. For each possible charge of 
the quasiparticles, m^^, with m = 0,l,...,/c, we only 
consider the particles with the lowest scaling dimension. 

These are given by c^^^e^^'^/^V^^^). There are two 
contributions to the scaling dimension: the parafermion 
field contributes A^rn = and the charge sec- 

tor contributes A0(m) = ^^j^2) •> giving ^ total scaling 
dimension Aj = (^^2)' which is therefore proportional 
to the charge of the excitations. 

We would like to create a situation in which we have 
one non-trivial quasiparticle that is degenerate with the 
vacuum, and an appreciable gap to the other types of 
quasiparticle excitations. The first thing we could try to 
do is to lower the energy of the charge e/{k + 2) fun- 
damental quasiparticle by means of an added potential, 
such that it becomes degenerate with the vacuum. How- 



ever, we just saw that adding such a potential will actu- 
ally create a set of + 1 degenerate states. To circumvent 
this problem, we will assume that, in addition to such a 
potential, there is also a charging energy proportional to 
q^^ where q is the charge of the excitation. Effectively, 
this modifies the amplitude of the quadratic contribution 
to the scaling dimension, coming from the charge part^^. 
Thus, by adding the charging energy, and the energy as- 
sociated with a suitable potential, we indeed can create 
the situation of two degenerate lowest lying states (one 
being the vacuum), separated from the others by a gap. 

We have just argued that we can consider a chain of 
itinerant anyons, consisting of vacancies with quantum 
numbers (/, 0) and fundamental quasiaparticle excita- 
tions with quantum numbers {^\^e/{k + 2)). Under fu- 
sion, the electric charge is merely additive, and we will 
therefore concentrate on the non-Abelian sector only. In 
the original anyonic chain models, the constituent anyons 
belong to the pure 811(2)^ theory. The anyonic systems 
we study can be maped to these by noting that the $j 
parafermionic field carries spin j = 1/2 SU(2)/c topo- 
logical charge, together with some Abelian topological 
charges. 
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The Ising TQFT has the same fusion algebra as SU(2)2, 
but the a and j = | anyons have different scaling dimen- 
sions (ha = 1/16 and hi = 3/16, respectively. The scal- 
ing dimensions of the and j — 1 anyons are identical, 
h^ = hi = 1/2). 

For k odd, the restriction to the integer- valued topologi- 
cal charge (generalized angular momentum) can be made 
by using the map obtained by fusing with the topological 
charge i = f , namely j x | — f ~ i 5 which for /c = 3 maps 
I o 1 and 1^0. For more details, we refer to Refs. 72,75. 
For the bosonic states, the particle identified with the vac- 
uum is a boson. In the fermionic case, one cannot sim- 
ply identify the electron with the vacuum, because it is a 
fermion, which obviously has different braiding statistics 
than the vacuum (which is a boson). Considering fusion 
and braiding, one could instead simply identify pairs of 
electrons with vacuum. However, the resulting theory will 
not be modular, meaning the S-matrix is degenerate. This 
poses a problem when one wishes to define the theory on 
arbitrary surfaces, including the torus. A solution is to put 
each charge into a Z2 doublet, e.g. the vacuum and elec- 
tron form the vacuum doublet, and every charge together 
with the charge obtained by fusion with an electron form 
a doublet. Then the S-matrix of doublets is modular. In 
practice, one can take the fusion rules assuming identifica- 
tion of electron with vacuum. 



